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' Abstract 

QhI We discuss the definition and the energy evolution of scattering amphtudes with C- 



odd ("odderon") quantum numbers within the effective theory for the Color Glass 
Condensate (CGC) endowed with the functional, JIMWLK, evolution equation. 
We explicitly construct gauge-invariant amplitudes describing multiple odderon ex- 
^ ' changes in the scattering between the CGC and two types of projectiles: a color- 

■ singlet quark-antiquark pair (or 'color dipole') and a system of three quarks in a 

colorless state. We deduce the energy evolution of these amplitudes from the gen- 
eral JIMWLK equation, which for this purpose is recast in a more synthetic form, 
which is manifestly infrared finite. For the dipole odderon, we confirm and extend 
the non-linear evolution equations recently proposed by Kovchegov, Szymanowski 
and Wallon, which couple the evolution of the odderon to that of the pomeron, and 
predict the rapid suppression of the odderon exchanges in the saturation regime at 
high energy. For the 3-quark system, we focus on the linear regime at relatively low 
energy, where our general equations are shown to reduce to the Bartels-Kwiecinski- 
Praszalowicz equation. Our gauge-invariant amplitudes, and the associated evo- 
lution equations, stay explicitly outside the Mobius representation, which is the 
Hilbert space where the BFKL Hamiltonian exhibits holomorphic separability. 
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1 Introduction 



Since the advent of the Bahtsky-Fadin-Kuraev-Lipatov (BFKL) equation 
[1,2] in the mid seventies, there has been significant progress in our compre- 
hension of high-energy QCD, and several theoretical approaches have been 
proposed which aim at a resummation of the energy-enhanced radiative cor- 
rections to high-energy processes in perturbative QCD. The BFKL equation 
is a leading logarithmic approximation (LLA), which allows one to resum to all 
orders corrections of the form (ug In s)" to the scattering between two colorless 
objects via the exchange of two gluons in the t-channel. As a result of this 
resummation, the bare two-gluon exchange is replaced by the BFKL pomeron 
(the sum of an infinite series of ladder diagrams of ordinary perturbation the- 
ory), or, equivalently, by two reggeized gluons which interact with each other. 
All the subsequent theoretical approaches proposed within perturbative QCD 
encompass the BFKL equation, and can be viewed as extensions of the latter 
towards increasing the complexity of the objects exchanged in the ^-channel, 
and also towards enlarging the limits of the LLA. 

The simplest object beyond the BFKL pomeron within perturbative QCD is 
the exchange of three interacting (reggeized) gluons in a symmetric color state. 
This object, which is negative (or "C-odd") under charge conjugation (C = 
— 1), represents the lowest order perturbative contribution to the odderon, the 
C~odd exchange which dominates the difference between the hadronic cross 
sections for direct and crossed channel processes at very high energies [3]. 
The evohition of the three-gluon odderon exchange with increasing energy in 
the LLA is described by the BKP equation, established by Bartels [4] and 
Kwiecinski and Praszalowicz [5], which amounts to a pairwisc iteration of 
the BFKL kernel (see also [6]). This equation can be immediately extended 
to describe the exchange of an arbitrary number n > 3 of reggeized gluons 
with pairwise BFKL interactions [7,4,8,9,10]. The resulting formahsm, also 
known as the generalized leading logarithmic approximation (GLLA), resums 
all radiative corrections that involve the maximally possible number of energy 
logarithms Ins for a given number of exchanged gluons. At the moment, two 
exact solutions of the BKP equation for odderon evolution are available [11,12], 
and the subject continues to be under intensive debate [13,14] (see also the 
recent review paper [15] and the discussion below). 

In the formalisms described so far, the number of gluons in the t-channel 
remains fixed in the course of the evolution. This is probably a good ap- 
proximation in some intermediate kinematical region, but it fails to describe 
two interesting physical situations: First, it does not incorporate correctly the 
fluctuations in the number of gluons, as resulting from processes in which one 
(reggeized) gluon splits into two, or, more generally, a n-gluon state evolves 
into a (n -|- m)-one, with m > 1. Such processes are especially important 
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in the dilute regime at relatively large transverse momenta (for a given en- 
ergy), where gluon splitting is the main process through which higher-point 
correlations get built [16]. Second, the approximation in which the number 
of i-channel gluons is fixed cannot describe recombination processes in which 
(reggeized) gluons merge with each other, thus reducing the gluon density. 
Such processes are important in the high-energy regime where the gluon den- 
sity becomes large enough (due to BFKL evolution and to the splitting pro- 
cesses alluded to above) to enhance recombination processes, which are then 
expected to lead to gluon saturation [17,18,19]. The inclusion of saturation is 
also necessary, for consistency, in studies of the unitarization of the scattering 
amplitudes, except for some exceptional kinematical configurations [20]. 

The simplest approach including gluon splitting in the framework of BFKL 
evolution is the color dipole picture developed by Mueller [20,21]. This picture 
is valid at large Nc, and describes pomeron multiplication via vertices at which 
one (BFKL) pomeron splits into two. A more ambitious program, which is not 
restricted to the large-A^c approximation, is the extended generalized leading 
logarithmic approximation (EGLLA), initiated by Bartcls [22], in which the 
gluon number changing vertices are explicitly computed in perturbative QCD 
(see Refs. [23,24,25,26,27,28,29] for further developments along this line and 
Ref. [15] for a review). By using such vertices, evolution equations allowing for 
gluon sphtting have been written down in Refs. [24,27,29]. Also, the equiva- 
lence between the triple pomeron vertex in the dipole picture [21,30] and the 
one generated by EGLLA at large Nc [23,25] has been verified in Refs. [24,28]. 

So far, the only formalism allowing for the systematic inclusion of gluon merg- 
ing in the high-energy evolution is the Color Glass Condensate (CGC) [31], in 
which the reggeized gluons are replaced by classical color fields whose correla- 
tions get built in the course of the evolution. But the corresponding evolution 
is non-linear : the new gluons radiated at one step in the evolution (the analog 
of the 'rungs' in the BFKL ladders) are allowed to scatter off the classical color 
fields generated in the previous steps, and this is the mechanism leading to 
gluon merging. Because of the non-linear effects, the evolution couples n-point 
functions with various values of n, and can be most compactly summarized 
as a functional Fokker-Planck equation for the weight function describing the 
correlations: the Jalilian-Marian-lancu-McLerran-Weigert-Leonidov-Kovner 
(JIMWLK) equation [32,33,34]. Alternatively, and equivalently [35], the evolu- 
tion can be formulated as an hierarchy of equations for scattering amplitudes 
— the Balitsky equations [36] — , in which unitarity is manifest. Note how- 
ever that gluon splittings are not included in the JIMWLK equation [16] ; this 
is obvious from the fact that, in the dilute, or weak-fi.eld, limit, this equa- 
tion reduces to an evolution in which the number of gluons in the t-channel 
stays constant [33,37]. An extension of the JIMWLK-Balitsky evolution which 
includes pomeron splitting has been proposed only very recently [16,38]. 
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As it should be clear from this succinct presentation, the various formalisms 
proposed so far in perturbativc QCD at high energies are quite different from 
each other, and the correspondences between them are not always transpar- 
ent. We know for instance that all these approaches reproduce the Balitsky- 
Kovchegov (BK) equation [36,39], which is the simplest non-linear generahza- 
tion of the BFKL equation, but only in the sense of a mean field approximation 
that has been recently challenged [16,40,41,42]. But the relation between the 
correlations (i.e., the n-point functions with n > 2) generated by the differ- 
ent approaches is much less understood. For instance, it has been shown only 
recently, by Kovchegov, Szymanowski and Wallon [43], that the perturbative 
odderon can be accommodated within Mueller's dipole picture [20], and that 
the corresponding solution coincides with the Bartels-Lipatov-Vacca (BLV) 
solution [12] to the BKP equation. 

In particular, in the regime where saturation effects can be neglected, one 
expects the CGC formalism and the more traditional approaches like GLLA 
to be equivalent with each other, but this has never been verified beyond the 
example of the 2-point function (i.e., of the BFKL equation). With this paper, 
we would like to make one more step towards elucidating this correspondence, 
by establishing the equivalence between the two approaches at the level of odd- 
eron exchanges (i.e., for a 3-point function). Specifically, we shall demonstrate 
that, in the weak-field limit, the JIMWLK evolution of the C-odd three-gluon 
exchanges reduces to the BKP equation, as expected. 

But recovering the BKP equation from the CGC formalism is not the main 
purpose of the present analysis, but only a pretext for it. The CGC is the 

theoretical framework par excellence for a study of high-energy scattering 
and evolution in QCD near the unitarity limit, yet the odderon problem has 
never been addressed in this formalism before. Thus, a substantial fraction 
of the subsequent analysis will be devoted to the proper formulation of the 
odderon exchanges in the framework of the CGC, and to the derivation of the 
corresponding evolution equations from the general, JIMWLK, equation. This 
study of the odderon should be a good starting point towards understanding 
the multi-reggeon dynamics within the CGC formalism. 

Our study will also emphasize some essential differences between the CGC for- 
malism and the perturbative approach based on the BFKL Hamiltonian: The 
latter is adapted to the description of a single scattering via the exchange of a 
composite object — pomeron, odderon, or, in general, a system of n reggeized 
gluons — which evolves with increasing energy. It relies on "/cr-factorization" 
(see, e.g., [44]) to separate the dynamics in the transverse plane from that in 
the longitudinal direction, and express a scattering amplitude as the convolu- 
tion of an universal Green's function, which describes the exchanged object, 
with the process-dependent impact factors, which connect this object to the 
external particles. From the above, one sees that the calculation is most nat- 
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urally carried on in momentum space. 



By contrast, in the CGC formalism — which is specially tailored to describe 
unitarity corrections — , single and multiple scatterings are treated on the 
same footing, namely they are resummed in process-dependent, and gauge- 
invariant, scattering amplitudes, which arc computed in the eikonal approxi- 
mation, and thus are naturally constructed in coordinate space. There is no 
/cy-factorization any longer, nor universal Green's functions: the longitudinal 
and transverse dynamics are tied up together in Wilson lines, which describe 
the eikonal scattering of the elementary particles which compose the projectile 
(the external object which scatters off the CGC, identified as the target). 

These differences explain some of the subtleties that we shall meet when trying 
to compare results for the odderon in the two approaches. On one hand, the 
odderon is described by the universal Green's function of three reggeized glu- 
ons, which obeys BKP equation in momentum space. On the other hand, the 
CGC scattering amplitudes depend upon the specific process at hand (they 
include the impact factor of the projectile) and obey non-linear evolution 
equations written in coordinate space. (In general, these are not closed equa- 
tions, but just a part of Balitsky's hierarchy [36].) Still, in the weak-field, or 
single-scattering, approximation, in which the evolution equations become lin- 
ear, they must contain the same non-trivial information as the BKP equation, 
whatever is the process under consideration. 

The authors of Ref. [43] have met with a similar difficulty when trying to 
compare the C-odd scattering amplitude of a dipole with the standard BKP 
odderon. In that case, they have been able to do so by using the respective 
solutions, which arc explicitly known. Here, we shall follow a more general 
strategy, which applies to arbitrary processes, including those where the evo- 
lution equations are too complicated to be solved exactly. Namely, by inspec- 
tion of two specific processes, wc shall be able to identify the analog of the 
universal odderon Green's function in the weak-field limit of the CGC for- 
malism, and show that, when properly defined, this quantity obeys indeed the 
(coordinate version of the) BKP equation. As we shall momentarily explain, 
this CGC approach to the BKP equation not only establishes a correspon- 
dence between the two formalisms, but also reveals some new insights about 
the BKP equation itself. 

The two specific processes that we shall consider are the CGC scattering with 
a quark-antiquark color dipole (a sub-process of the virtual photon — CGC 
scattering) and that with a colorless 3-quark system (a simple model for a 
baryon). For both cases we start by constructing the general, non-linear, am- 
plitudes which describe multiple odderon exchanges (these turn out to be the 
imaginary parts of the respective ^-matrix elements, themselves expressed 
in terms of Wilson lines), and then expand these amplitudes in the limit 
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where the CGC field is weak. After this expansion, both amphtudes reduce to 
(gauge-invariant) Unear combinations involving a three-gluon Green's func- 
tion in a totally symmetric color state. Clearly, this Green's function is a 
natural candidate for the BKP odderon in the CGC formalism. This inter- 
pretation is, however, hindered by the fact that the CGC Green's functions 
are gauge-variant objects, for which the JIMWLK equation predicts infrared 
singularities (to be contrasted with the BKP equation, which is infrared safe). 
Although physically harmless — as they cancel in the gauge-invariant am- 
phtudes — , these singularities comphcate the correspondence with the BKP 
approach. 

At this point comes one of the main new technical developments in this paper: 
We show that the JIMWLK Hamiltonian [33,34] can be rewritten in a new 
form, which is manifestly infrared finite (the original kernel in the transverse 
space is replaced by the dipole kernel [20], which decays much faster at large 
distances). When acting on gauge-invariant quantities, this new Hamiltonian 
is equivalent with (but simpler to use than) the original one, in the sense of 
generating the same evolution equations. But the new Hamiltonian generates 
infrared-finite equations also for the gauge-variant Green's functions, and thus 
allows us to introduce the latter in a mathematically well-defined way. With 
this prescription, the equation satisfied by the CGC odderon Green's function 
turns out to be the same as the Fourier transform to coordinate space of the 
BKP equation, as we shall check explicitly. 

But this Fourier transform reserves some more surprises, as it could be an- 
ticipated from the fact that our equation in coordinate space is not exactly 
the same as the coordinate-space version of the BKP equation that is usually 
written in the literature^ (see, e.g., [15]). Rather, the two equations coincide 
with each other only if we require our CGC Green's function, which in gen- 
eral is a totally symmetric function of three transverse coordinates, to vanish 
whenever two coordinates become identical. This property is sometimes re- 
ferred to as "the Mobius representation" (see, e.g., [28]), and is interesting 
in that, when restricted to functions having this property, the BFKL Hamil- 
tonian is conformally invariant [2] and exhibits holomorphic separability [8]. 
This mathematical simplification has led [9,10] to a powerful analogy between 
the BKP odderon problem (and, more generally, the problem of multi-reggeon 
exchanges in the limit of a large number of colors) and an integrable Heisen- 
berg spin chain. In particular, the first exact solution to the BKP equation 
has been found, by Janik and Wosiek [11], by exploiting this analogy. 



^ At the technical level, the difference originates in some ambiguities in the form 
of delta-functions which appear when Fourier transforming the momentum-space 
BKP equation to coordinate space, and which are generally interpreted in the sense 
of the Mobius representation [28]; that is, these delta-functions are simply ignored. 
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But although natural for the pomeron exchange (i.e., at the level of the 2-point 
function), where it entails no loss of generality, the restriction to the Mobius 
representation is not so natural for the higher n-point functions (n > 3), as 
recently emphasized in Ref. [28]. For instance, the other known solution to 
the BKP equation, due to Bartels, Lipatov, and Vacca [12] , which dominates 
at high energy and is perhaps of more relevance for the phenomenology (as it 
couples to a virtual photon), lies outside the Mobius representation. 

Similarly, the use of the Mobius representation does not appear to be natural 
in the CGC formalism either. In fact, our both examples of gauge-invariant 
scattering amplitudes lie outside this representation: For the dipole case, there 
is no coupling to this functional space, as well known [12,43], whereas for the 
3-quark case, this property is excluded by the initial conditions. We conclude 
that, at least for the problems that we shall discuss, the BKP equation must 
be solved in a Hilbert space more general than the Mobius representation. 

But our analysis below will not be confined to the weak-field limit and its re- 
lation with the perturbative QCD approaches. As repeatedly emphasized, the 
CGC is a formalism for multiple scattering, in which non-linear amplitudes 
and the corresponding evolution equations are straightforward to construct. 
As an illustration, we shall derive the general evolution equations for the 
scattering amplitudes describing C-odd and, respectively, C-even exchanges 
in the dipole-CGC scattering (the non-linear generalizations of the odderon 
and, respectively, pomeron exchanges). In fact, these equations will be ob- 
tained by simply separating the real part and the imaginary part of the first 
equation in the Balitsky hierarchy [36]. Interestingly, the non-linear terms in 
these equations are found to couple the odderon and pomeron evolutions. In 
the mean field approximation in which the non-linear terms are assumed to 
factorize, the equation for the C-odd ampUtude reduces to a non-linear equa- 
tion originally proposed in Ref. [43] . Our analysis of this equation will confirm 
the conclusion [43] that the odderon exchanges are strongly suppressed by the 
unitarity corrections, and will allow us to deduce the mathematical law for 
this suppression. For the 3-quark system, we shall not write down the corre- 
sponding non-linear equation (since this appears to be too complicated to be 
illuminating). Rather, we shall rely on the relation between the corresponding 
C-odd amplitude and the respective one for the dipole to conclude that, in the 
weak-field regime, the dominant increase with the energy should be controlled 
again by the BLV solution [12], so like for the dipole case [43]. 

The plan of the paper is as follows. In Sect. 2 we give a general proof that the 
JIMWLK evolution of gauge-invariant observables is free of infrared problems, 
and we deduce an alternative form of the JIMWLK Hamiltonian which makes 
infrared finiteness manifest. In Sect. 3, we consider the weak-field limit of the 
JIMWLK evolution, and show that the use of the new Hamiltonian allows 
one to introduce well-defined CGC Green's functions. In Sect. 4 we construct 
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the general amplitudes describing niTiltiplc C-odd exchanges for a color dipole 
and a 3-quark system. Then, in Sects. 5 and 6, we deduce the corresponding 
evolution equations, after having introduced first the odderon Green's function 
in the CGC. Finally, in Sect. 7 we discuss the connection to the BKP equation. 



2 The JIMWLK equation with the dipole kernel 

In this section, we show that the JIMWLK equation can be considerably sim- 
plified when its action is restricted to gauge-invariant correlation functions, 
such as the scattering amplitudes. The resulting equation is still a functional 
differential equation, but with a different kernel in transverse space — the 
dipole kernel — , which has a more rapid fall-off at large distances, and there- 
fore makes it easier to check that the evolution is free of infrared singularities. 

2.1 The JIMWLK equation 

The CGC formalism [31] is an effective theory for the small-x gluons in the 
light-cone wavefunction of an energetic hadron. In this formalism, the gluons 
with small longitudinal momenta, or small values of x, are described as the 
classical color field radiated by 'color sources' (gluons and valence quarks) with 
higher values of x, which are 'frozen' by Lorentz time dilation in some random 
configuration. Accordingly, the color fields at small-x are themselves random, 
with a distribution specified by the 'weight function' VFr[<^] (a functional prob- 
ability density). Here, r = ln(l/a;) is the rapidity, and a = Q;''(a;^,cc) is the 
light-cone component of the color gauge field, and is the only non-trivial com- 
ponent in a suitable gauge (the 'covariant gauge'; see below). Note that this 
field depends upon the light-cone longitudinal coordinate ^ = (t — z)/ v^, 
and upon the transverse coordinates x = {x,y), but not upon the (light-cone) 
time x~^ = [t + z)/ \/2, in agreement with the 'freezing' property mentioned 
above. 

All the interesting physical quantities are expressed as operators built with a, 
say and the corresponding expectation values are obtained after aver- 

aging over the random field a: 



Whereas, by itself, the weight function is a non-perturbative object, its 

evolution with decreasing x (or increasing energy) can be computed in pertur- 

^ We assume that the color glass moves in the positive z direction. 
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bation theory, at least in the high energy regime where the intrinsic 'saturation 
momentum' Qs{x) (which increases as a power of 1/x [17,58]) is hard. The 
corresponding evolution has been computed in the non-linear generalization 
of the leading logarithmic approximation [32,33], which allows one to extend 
the BFKL resummation [1] in the high density region at saturation. In this 
resummation, the radiative corrections enhanced by the logarithm of the en- 
ergy Ins ~ lnl/,T, and the non-linear effects involving the classical field a, 
are all treated on the same footing, as effects of order one. The result of this 
calculation [33] is a second-order, functional, differential equation for WV[a], 
which is known as the JIMWLK equation [32,33,34] and reads: 



^WAa] = HWAa] ^ \ j y)j^^M, (2-2) 

where the subscript xy on the integral sign denotes the integration over the 
transverse coordinates x and y. The kernel r]"-''{x, y) is a functional of a, upon 
which it depends via the Wilson hnes V{x) and V\x) built with a = a"-T"- 
in the adjoint representation: 

= (i-Kiv;)''^(i-v;tv^,)'\ (2.3) 

with the following transverse kernel: 



]C{x,y,z) 
and, e.g.. 



{x- z)-{y- z) 
(x - zf{z - yf 



(2.4) 



Vl = 0(a;) = Pexp (^ig j dx-a^ix', a;)T"^ , (2.5) 

where P denotes path-ordering in x~ , and the integration over runs over the 
longitudinal extent of the hadron, which increases with r : When decreasing 
X, we include in the effective theory gluon modes with smaller longitudinal 
momenta, which by the uncertainty principle are localized at larger values of 
x^ . The functional derivatives in Eq. (2.2) act on the color field created in 
the last step of the evolution (i.e., in the rapidity bin (r, r + dr)), which is 
therefore located at the largest value of x~ . Thus, the action of the derivatives 
on Wilson hnes like (2.5) reads as follows: 

^^.,^«(.-,)r-v-., J|^ = -.,,<^)(.-,)v;t". (2.6) 
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By taking a r-dcrivative in Eq. (2.1) and using the JIMWLK equation (2.2), 
one can easily deduce the following evolution equation for a generic observable: 



-7- 

For O to represent a physical observable, this must be a gauge-invariant op- 
erator, or, more precisely, the expression of such an operator when evaluated 
in the covariant gauge in which A^^^x) — 5'*+Q;a(x~, x). In such a case, one has 
seen on specific examples that the evolution described by Eq. (2.7) is infrared 
safe [33], and in what follows we shall give a general proof in that sense. 

As a simple example, consider the scattering between the CGC and an external 
'color dipole' : a quark-antiquark pair in a colorless state. The corresponding 
(S-matrix operator can be computed in the eikonal approximation as: 



= —ti{VX)^S{x,y-a), (2.8) 

where is a Wilson line in the fundamental representation (as obtained by 
replacing T" — > t"- in Eq. (2.5)), and represents the phase factor picked up by 
the quark while crossing the background field of the target. Similarly, Vy is 
the corresponding phase factor for the antiquark. Plugging this operator into 
Eq. (2.7), one obtains [33] the following equation (with ots = agNc/n): 



§^{i<VX))^ = ^J d'zM{x, y, z) [j^i^V^V^MVlVy) - tr(V;tV,))^ , 

(2.9) 

after a rather lengthy calculation in which many terms which appear at inter- 
mediate steps cancel with each other. This is not a closed equation, but only 
the first equation in an infinite hierarchy originally derived by Balitsky [36]. 
The kernel appearing in this integral has been generated as: 



M{x,y,z) = _ ^ iC^xz + iCyy;^ - 21C^y^ , (2.10) 

[X z) [z y) 

and is recognized as the dipole kernel [20]. Note that the poles in this kernel 
at 2; = a; and z = y are actually harmless because the operator within the 
brackets vanishes at these points. In fact, it is easy to check on Eq. (2.3) 
that such 'ultraviolet' (i.e., short-distance) poles cancel already in the general 
evolution equation (2.7), irrespective of the nature of the operator O. 
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The crucial new feature of the dipole kernel (2.10) as compared to the original 
kernel (2.4) in the JIMWLK equation is that this new kernel falls off much 
faster at large distances: /A{x,y,z) ^ when x,y, which is enough 
to ensure the convergence of the integral in the r.h.s. of Eq. (2.9). That is, this 
equation and, similarly, all the higher equations in the Balitsky hierarchy [36], 
are infrared safe. As we shall explain in the next subsection, this is related to 
the fact that the corresponding operators are gauge invariant. 

2.2 The dipole JIMWLK equation and the finiteness conditions 

A brief inspection of the original JIMWLK equation, cf. Eqs. (2.2)-(2.4), re- 
veals that, for a generic operator O, there is a priori no guarantee that the 
corresponding evolution equation (2.7) should be infrared safe. Indeed, at large 
distances z x, y, the transverse kernel (2.4) decays only like }Cxyz ~ 
so Eq. (2.7) may develop logarithmic divergences (~ / d^z/ for z oo). In 
the calculation leading to (2.9), such divergences have actually canceled among 
various terms. In what follows, we shall argue that such a cancellation will al- 
ways hold for the gauge-invariant operators. In the present context, where 
a gauge-fixing condition has been already chosen, what we mean by that is 
the invariance under the residual gauge transformations. An interesting con- 
sequence of this discussion will be that, for such operators, the corresponding 
evolution equation (2.7) can equivalently be rewritten in a simpler, and man- 
ifestly infrared finite form, which involves the dipole kernel (2.10): 

^<^>-<«*^>--i^/(J^f)i(f^ (2.11) 

Our argument will be constructed as follows: First, we shall use Eqs. (2.3) and 
(2.7) to deduce the "finiteness conditions" that some operator O must satisfy 
in order for the corresponding evolution equation to be infrared safe. Then, 
we shall use these conditions to rewrite Eq. (2.7) in the form of Eq. (2.11). 
Finally, we shall demonstrate that the "finiteness conditions" are equivalent 
to invariance under the residual gauge transformations in the covariant gauge. 

The conditions for infrared finiteness are easily written down once one realizes 

that the dangerous terms are those terms in Eq. (2.3) which do not involve 
either V^^or^Vj". Indeed, after rewriting (1 - t/jt4)(l - V}.Vy) = 1 + V^Vy - 
V^Vz — V^Vy, it becomes clear that the contribution of the first two terms, 
1 -I- Vj.Vy, to Eq. (2.7) is potentially divergent: 
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xyz 



{x-z).{y-z) 5 , 5 

{x - zf{z - yf Sa^ixy ^ ^ y' Sa^^iy) 

f T4^-^(^ + ^l^yT'T4^-^^^ (forz^oo), 
J z^ J Sa":(x)^ ^ofiAy) 

whereas the contributions of the other terms should be finite, because corre- 
lations involving Vz or are expected to vanish when, e.g., \z — x\ oo. 
Clearly, for the divergence in Eq. (2.12) to go away, the coefficient of the di- 
vergent z-integral there must vanish, which in turn implies the two following 
finiteness conditions: 



I 
J 



It is easy to check that the following 2n-point operators constructed from the 
Wilson lines in the fundamental representation 

On = tr{MiM2---Mn), M, = Vl Vy^ , (2. 14) 

(and arbitrary linear combinations and products of them 0^ ■ ■ ■) satisfy 
the finiteness conditions Eqs. (2.12) and (2.13). For example, consider the 
simplest case n — 1, i.e., the dipole operator (2.8). The first condition (2.12) 
is trivially satisfied because 



j^fiVX) = ^9 {S^'\v -X)- S^'\v - y)} tr{VXn 

which indeed vanishes after integration over v, while the second condition 
(2.13) can be verified with the help of the formulae {V)f = 2tr(rK^VJ') and 
tr(iM)tr(rS) = |tr(AS) - ^tr{A)tT{B). 

Let us now establish the dipole form of the JIMWLK equation, Eq. (2.11). We 
first decompose the original kernel (2.4) so that the dipole kernel is separated 
out (cf. Eq. (2.10)): 



{x- z)-{y- z) ^ 1 
{x-z)\y-zy~ 2 



{x-zY{y-zf {x-zy {y-zf 



Note that this is in fact the separation of infrared divergences: The last two 
terms will generate logarithmic divergences at large z when the JIMWLK 
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Hamiltonian is applied to a generic operator. However, it is easy to show that 
if the operator satisfies the finiteness conditions, the contributions coming 
from the last two terms \/{x — zY and l/{y — zY vanish. Thus, for such 
operators, we can replace JC{x,y,z) — > —{l/2)M.{x,y,z) in Eq. (2.3), and 
at the same time move the functional derivative S/6a^{x) to the right of the 
Wilson lines. Indeed, the relevant commutator, namely. 



vanishes when multiplied by the factor [x — yY in the dipole kernel. This 
establishes Eq. (2.11). 

Besides being conceptually more appealing (as infrared finiteness is now man- 
ifest), the rewriting of the JIMWLK equation as in Eq. (2.11) considerably 
simplifies the manipulations leading from Eq. (2.11) to final evolution equa- 
tions like Eq. (2.9). The dipole kernel, that we expect in the final equations, is 
already present there, and many of the terms generated at intermediate steps 
when the calculations are based on the original equations, Eqs. (2.7)-(2.3), are 
simply absent when the calculations start with Eq. (2.11). For the more com- 
plicated operators that we shall encounter in the next sections, the usefulness 
of the dipole JIMWLK equation is indeed appreciable. 

2.3 Physical meaning of the finiteness conditions 

The operator (2.14) which has been seen to obey an infrared-finite evolution 
equation is the covariant-gauge expression of a gauge invariant operator. In- 
deed, one can rewrite this operator as the trace of a closed Wilson loop, for 
which gauge symmetry is manifest. To that aim, note that the end points at 
longitudinal infinity {x' — ±oo) of two adjacent Wilson fines (say, V^. and Vy^) 
can be connected by a Wilson line in the transverse direction, which is simply 
unity in the present gauge, in which = 0. Therefore, one can connect the 
end points of all the Wilson lines which enter the trace in Eq. (2.14) can be 
connected in such a way to form a single, closed, Wilson loop (see Figure 1). 

This observation suggests that there is an intimate relation between the in- 
frared finiteness and gauge invariance, that we would like to clarify now. To 
this end, we first notice that the differential operators which enter the finite- 
ness conditions, namely 

6.W^-£/rf^a=^. (2.16) 
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Fig. 1. A closed Wilson loop for ti(VjVyVjVuj) as an example of the operator (2.14). 
Solid lines are Wilson lines along the longitudinal directions which are manifest on 
the operator, and dashed lines are Wilson lines which connect two end points at 
longitudinal infinity x~ = ±00, but are in fact unity in the covariant gauge. 

^rK] ^ -u^^ J d'x{V.rjJ^^ , (2.17) 
are, respectively, the generators of left and right, global, color rotations: 



V^^n^Vl (2.18) 
V^^^V^nl (2.19) 

where 1^l/r = ^'^p{w^l/r^'^) ^ constant SU(A''c) matrix. Indeed, by using the 
formulae analogous to Eq. (2.6), one can explicitly check that the infinitesimal 
color rotations are given by 



SlV^ = QlIoolM = WuJlVI (2.20) 
SrV^ = GkMV^ = -tgV^u;^. (2.21) 

Thus, for any operator O biiilt with Wilson lines, the finiteness conditions 
(2.12) and (2.13) are tantamount to the conditions of invariance under such 
global color rotations. It is trivial to check that the operator (2.14) is indeed 
invariant under the color rotations (2.18) and (2.19). 

Now, with their one-sided action, the color rotations (2.18) and (2.19) look a 
priori different from the ordinary gauge transformations of the Wilson lines. 
Still, as we explain now, they are in fact the residual gauge transformations 
with respect to the 'covariant gauge' in which the CGC theory is usually 
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formulated (see [31] for details). We recall here that, in this gauge, the vector 
potential has only one non-zero component, the light-cone field = a^, 
which is independent of x^. This structure of the field is preserved by a gauge 
transformation A^^ — > Q{A^ + in which the gauge function Q depends 

only upon x". Thus a residual gauge transformation amounts to 



A+ ^ n{x-) (^A+ + ^d+j n\x-), nix-) = e'^'^"(^")*", (2.22) 
which induces the following transformation for the Wilson fine : 



n{x- = oo) V^n^ix- = -oo). (2.23) 

So far, the gauge function uj{x~) is arbitrary. The "global" color rotations 
introduced in Eqs. (2.18) and (2.19) are now obtained as the two (independent) 
special gauge transformations in which uj{x~) — > (and thus Q{x~) ^ 1) at 
either x~ — > — oo (for the left rotation (2.18)) or x~ — > +oo (for the right 
rotation (2.19)). This estabhshes the interpretation of the finiteness conditions 
as the conditions for gauge symmetry ^ . 

Finally, let us show the infinitesimal transformation for the gauge field inte- 
grated over the longitudinal direction: q;"(cc) = / dx~ a"" {x~ , x) . As we will 
discuss shortly, this is a natural variable in the weak-field regime. Since 
^4+ = a{x~, x) transforms as in Eq. (2.22), its infinitesimal change is given by 
5a{x~,x) — d'^u;{x~). Therefore, a"'{x) transforms as 



a^ix) — ^ a''{x)+ j dx-d+co^ix-) = a''{x)+C: (2-24) 

— oo 

with ^" = u!"'{x~ — oo) — a;"(x~ = — oo) being a pure number. This transfor- 
mation law will be useful in checking the gauge invariance of various operators 
in the weak-field limit. 



^ Note also that the left and right generators form two independent SU(A^c) algebra, 
as is expected for gauge transformations. More precisely, the differential operators 
J£(a;) = -ij^, and J^{x) = j^iVx)''^^^, which enter the definitions of the left 

and right generators, satisfy the following algebra [Sxy = S^"^^ {x — y)] 

[Jlix), Jliy)] = ir''Jl{x)6,y, [rnix), jUv)] = ir"'J%{x)S,y, 
and commute with each other [J£(a;), J^^iy)] = 0. 
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3 WeEik— field regime and Green's functions in the CGC 



In order to make contact with previous evolution equations proposed in per- 
turbative QCD at high energy, like the BFKL [1] and BKP [4,5] equations, 
which neglect saturation effects and therefore apply only in the dilute regime 
at relatively high transverse momenta (well above the saturation scale), it 
is convenient to consider also the weak-field limit of the JIMWLK equation 
(here, with dipolar kernel). When ga <^ 1, one can expand the Wilson lines 
in Eq. (2.3) to lowest nontrivial order, to obtain 



(1 - V^V^y^ ^ tg{a%x) - a^(z))(n^„, (3.1) 
where 

a^ix) = j dx- a''{x-,x) = < (3.2) 

is the effective color field in the transverse plane, as obtained after integrating 
over the longitudinal profile of the hadron. Then, the kernel 7] becomes simply 
quadratic in a J, and the evolution equation (2.11) reduces to 

9 ,,^\ 9"^ I' - Vf 



xyz 



{x - zf{z - yf 

For consistency with the above manipulations, the observable O — which 
in general is some operator built with Wilson lines (see, e.g., Eq. (2.14)) — 
must be correspondingly expanded in powers of ga. Note that the relevant 
order for this expansion depends upon the operator at hand, and also upon 
the channel that we consider for scattering (see Sect. 4). For instance, for 
the dipole operator shown in Eq. (2.8), the lowest nontrivial contribution is 
obtained after expanding the Wilson lines up to second order in ga : 

V^[a] ^1 + ig J dx-a^ix', x)^ (3.4) 

2 

~Y I I '^y^^"^^'^ x)a\y-, x) [e{x- - y')rt^ + e{y- - x-)tH"'\. 

Clearly, to this order, the ordering of the color matrices Q;"(a;^, a;)t" in x^ starts 
to play a role in the expansion of the Wilson lines. But this ordering is still 
irrelevant for the computation of the dipole (S-matrix to lowest order, because 
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of the symmetry of the color trace in Eq. (2.8) : tr{tn^) = | = tr(i^f"). 
Namely, one finds: 



S{x,y;a)^—tr{V^Vy) 1 




(3.5) 



where the final expression involves only the two-dimensional field a^, cf. 
Eq. (3.2). The latter property turns out to be shared by all the operators 
for weak (single) scattering that we shall discuss throughout this work. This 
feature, together with the structure of the simplified evolution Hamiltonian 
manifest on Eq. (3.3), enable us to consider the weak field approximation to 
the CGC effective theory as being restricted to the Hilbert space of the func- 
tions built with a^. Onto this space, the functional derivatives in Eq. (3.3) can 
be trivially replaced by 5/Sa%. Thus, in this approximation, the longitudinal 
structure of the color field becomes irrelevant. 

Note that the weak-field Hamiltonian in Eq. (3.3) is quadratic both in a 
and in the functional derivative with respect to a. Thus, when acting on the 
n-point function {a{xi)a{x2) ■ • • Oi{xn))T, this Hamiltonian does not change 
the number n of fields. This has an important consequence for the weak-field 
evolution described by Eq. (3.3): During this evolution, the number of gluons 
in the t-channel remains fixed, familiar to the 'multi-reggeons' approaches 
based on BFKL evolution [7,4,5,8,10], to which we shall eventually compare 
the present formahsm. In particular, this imphes that the JIMWLK evolution 
in the weak-field regime is linear. 

It is already known [32,33] that, in the weak field limit, the JIMWLK evolution 
of the dipole -S-matrix reduces to the corresponding BFKL equation [1] . This 
can be easily checked on the first Bahtsky equation, Eq. (2.9), by first rewriting 
this equation in terms of the dipole scattering amplitude, 



and then linearizing the ensuing equation with respect to T, which is formally^ 
appropriate in the weak scattering regime where |T| <^ 1. But for the following 
developments in this paper it is still instructive to give a rapid derivation of the 
BFKL equation, by using directly the weak field approximation in Eqs. (3.3) 
and (3.5). Specifically, in this approximation —i{T{x, y))r — {N{x, y))T, with: 



^ We ignore here the subtleties associated with fluctuations in the dilute regime 
which in general render such a linearization illegitimate even when KT),-! <C 1; see 
Ref. [16] for details. 



iT{x, y;a) = l- S{x, y;a) = l 



1 



(3.6) 
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with the following definition for the 2-point Green's function of the color fields 
in the dilute regime: 



Mx, y) = {a\x)a\y))r = U{y, x). (3.8) 

Note that, although colorless (it carries no open color indices), the object in 
Eq. (3.8) is still not gauge invariant: the residual gauge transformation for a^. 
consists in the constant shift — > a" + [cf. Eq. (2.24)], but fT-{x,y) is 
not invariant under this operation. On the other hand, the linear combination 
yielding the scattering amplitude (3.7) is clearly invariant, as it should (since 
Eq. (3.7) has been obtained after a consistent expansion of the gauge-invariant 
operator (3.6)). Thus, a priori, one is allowed to use the dipolar form of the 
evolution equation, Eq. (3.3), for the physical amplitude ((a" — a^)^)T, but 
not also for the Green's function fr{x,y). Still, since Eq. (3.3) is linear in O, 
it is clear that the equation obeyed by {{a% — oty)'^)r is correctly obtained by 
separately evolving each of the Green's functions in the second line of Eq. (3.7), 
and then summing up the corresponding results. This argument shows that, 
in fact, it is legitimate to use the dipolar evolution equation (3.3) even for 
quantities which by themselves are not gauge invariant, so like the Green's 
function (3.8), provided these quantities are eventually used as building blocks 
in the construction of gauge-invariant observables. In such a case, the use 
of the dipolar Hamiltonian if^p can be seen as a convenient prescription to 
regulate the infrared singularities which would appear at intermediate steps 
when using the original JIMWLK Hamiltonian in Eqs. (2.2)-(2.3). 

Specifically, by using Eq. (3.3) for O — a'^oiy, one easily finds 



= (3.9) 

^{fAx^z) + U{y,z) - fr{x,y) - fr{z,z)). 

This equation is well defined both in the infrared (because of the rapid decay 
of the dipole kernel at large values of z), and in the ultraviolet (the short- 
distance poles of the kernel at x = z and y = z are actually harmless because 
they have zero residue). By contrast, the equation which is obtained by acting 
on a^Uy with the (weak-field version of the) original JIMWLK Hamiltonian ^ 
contains terms which are potentially singular at large distances. 



This equation can be found, e.g., as Eq. (3.10) in Ref. [45]. 
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By using Eqs. (3.7) and (3.9), one finds the coordinate space (or 'dipolar') 
version of the BFKL equation for {N{x, y))T, as expected: 



d_ 



{Nix,y))r = ^ld'z 



{x - yf 



(3.10) 



{x - zY{y - zf 



X 



[{N{x, z))r + {N{z, y))r - {N{x, y))r) . 



In this last equation, ultraviolet finiteness is ensured by the vanishing of the 
scattering amphtude at equal points, {N{x,x)) = 0, a property consistent 
with Eq. (3.7) and which reflects 'color transparency'. 

In what follows, we shall take Eq. (3.9) (supplemented with an appropri- 
ate initial condition) as the definition of the 2-point Green's function in the 
CGC formahsm. This situation illustrates a general feature of the present 
approach, namely the fact that, because of infrared complications, the defi- 
nition of Green's functions meets with ambiguities which disappear only in 
the construction of gauge-invariant quantities. Of course, these Green's func- 
tions are only intermediate objects, which are strictly speaking not needed: 
it is always possible to construct directly the equation obeyed by the gauge- 
invariant quantity of interest, which is then free of any ambiguity. The reason 
why it is nevertheless convenient to introduce such Green's functions, it is 
because these are the CGC analogs of the multi-gluon exchanges considered 
in the more traditional approaches to high energy QCD [7,4,5,8,10], to which 
we would like to establish a connection in the forthcoming sections. 



4 Odderon operators in the CGC 



With this section, we start our discussion of the odderon within the CGC 
framework. To start with, we shall construct the operators describing multi- 
ple odderon exchanges in the scattering between the CGC and a relatively 
simple projectile, such as a color dipole, or three quarks in a colorless state. 
The dipole-CGC scattering can be viewed as a sub-process of the diffractive 
scattering of a virtual photon on some dense hadronic target (the CGC), a 
process in which odderon contributions are expected, e.g., in the production 
of C-even mesons like rjc (see, e.g., [12,46]). As for the 3-quark system, this 
may be viewed as a crude 'valence quark' model of the baryon. 
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4-1 The dipole-CGC scattering 



Consider first the simplest case, namely the high energy scattering of a qq 
dipole off the CGC, and let us briefly outline the construction of the cor- 
responding S'-matrix. As explained in Sect. 2.1, we need to first compute 
the S'-matrix for a fixed configuration of the classical field a, and then aver- 
age over the latter. For the first step, we can use the eikonal approximation: 
S{x, y; a) = (out|in), where the transverse positions of the quark (x) and the 
antiquark (y) are the same in the in-coming and the out-going states. One 
can write, schematically, {i — 1, • ■ ■ , Nc is the color index) 

I in) ~ io) , I out) ~ vsr'(^«)V'r'(2/) |0) , 

where the appropriate normalization is understood. The relation between the 
in-coming and the out-going fields is found by solving (the high-energy ver- 
sion of) the Dirac equation (9_ — iga°'t"')t/j = for a given gauge field con- 
figuration a. This implies that = (y^)ijipj^ with the Wilson line in the 
fundamental representation. The S'-matrix becomes: 

S{x,y;a) = (out|in) = i-(V;t)^^(V^,)^W' = i-tr(KtK,), 

(we have restored the appropriate normalization), in agreement with Eq. (2.8). 
The physical S'-matrix is finally obtained after averaging over the random 
classical color field, cf. Eq. (2.1), an operation which also introduces the de- 
pendence upon the energy (i.e., upon r): 

Sr(x,y) ^ ^MVlVy)),. (4.1) 

Since non-linear in a. the above formula describes in general multiple ex- 
changes, which can be cither even, or odd, under the operation of charge 
conjugation C. To single out C-even ('pomeron') or C-odd ('odderon') ex- 
changes, one needs to project Eq. (4.1) onto incoming and outgoing states 
with appropriate C-parities. Since the charge conjugation for fermions is de- 
fined by CipC"-^ = — i(-?/'7°7^)'^, and CipC'^ = (— ■i7°7^-?/')^, it is clear that the 
eigenstates of C in the dipole sector are given by {ip{x)ip{y) ±ip{y)ip{x)) |0), 
where +(— ) sign yields the C-even(odd) state. These structures are natural 
because the charge conjugation is essentially the exchange of a quark and an 
antiquark. 

Taking the C-odd dipole state as the in-coming state (this is the state selected 
by the virtual photon wavefunction) , and the C-even dipole state as the out- 
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going state (this would be selected by a rjc meson), one obtains the following 
C-odd contribution to the ^'-matrix: 



S°''''{x,y) = (out, even|in, odd) = ^ (tr(yX) - tr(yjK))^ . (4.2) 

This allows us to identify the operator for C-odd exchanges in the dipole-CGC 
scattering ("the dipole odderon operator") as 



0{x, y) = i^HVlVy - V^V^) = - 0{y, x), (4.3) 

where the factor of i is introduced in order for this quantity to be real: indeed, 
since V and are unitary matrices, we have [tr(V^Vy)]* = txiV^Vx)- 

One can directly check, by using the transformation property of the gauge 
fields under charge conjugation. 



C A,C-' ^ -{A,f , (4.4) 

that the operator (4.3) is indeed C-odd: For a generic Wilson line V con- 
structed with A^, Eq. (4.4) implies 



CVC-^ = {V^f = V\ (4.5) 

so that CiT{V^Vy) C'^ = tr(yjT4), or, finally, C 0{x, y) C'^ = -0{x, y). 

Note that the C-odd contribution (4.2) is the imaginary part of the ^'-matrix 
element, or, equivalently, the real part of the scattering amplitude T, with 

{Oix,y)), = QmSrix,y), (4.6) 

which was to be expected. Correspondingly, the C-even, Pomeron exchange, 
amplitude, that we shall denote as N{x, y), is identified with the real part of 
the jS'-matrix: 



N{x, y)^l- ^t^iVX + ViVx), (4.7) 
{N{x,y))r^l-LSrix,y). (4.8) 

Operatorially, S = 1 + iT = 1 — N + iO. Note the obvious boundary conditions 
which follow from ^'^-(a;, x) = 1 (or directly from the definitions (4.3), (4.7)): 
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N{x,x) = 0{x,x) = 0. 



(4.9) 



Clearly, these conditions remain true after averaging over the random field a. 

From perturbative QCD, we expect the lowest order contribution to the odd- 
eron exchange to be represented by three gluons tied up together with the 
^abc gyiYihol, where d"''"^ = 2tr({t", t^jt*^) is the totally symmetric tensor. The 
operator d"^'^yl^(cc)A[^(?/)y4^(z) is indeed C-odd, as obvious from Eq. (4.4). 
Let us check that a similar structure emerges also from the CGC operator 
(4.3) when this is evaluated in the weak-field limit. The lowest non-trivial 
contribution to Eq. (4.3) is obtained by expanding the Wilson lines there up 
to cubic order in the field a in the exponent. (The terms quadratic in a cancel 
in the difference of traces in Eq. (4.3).) By collecting the remaining terms, one 
obtains: 



0{x, y) -^d^''^ {3(«a^ - alaial) + {a^a^al - a^^) j(4.10) 



As expected, this expression is cubic in a" with the color indices contracted 
symmetrically by the d-symbol. Note that, because of the symmetry properties 
of this symbol, the path-ordering of the Wilson lines in x~ has been irrelevant 
for computing 0{x,y). 

The hnear combination of trilinear field operators in Eq. (4.10) is gauge invari- 
ant by construction. To render this more explicit, let us rewrite this operator 
as 



This is manifestly invariant under a residual gauge transformation, which con- 
sists in a constant shift of the field: a"- a"" + (cf. Eq. (2.24)). In fact, an 
alternative way to construct the gauge-invariant linear combination appear- 
ing in Eq. (4.10) is to directly impose the (weak-field version of the) finiteness 
conditions (2.12) and (2.13) on the C-odd Green's function ^"^^(a^a^Q;^)^. 

4-2 The 3-quark-CGC scattering 

To describe the 3-quark colorless state, we shall use the following "baryonic" 
operator ^ e^^^ijj^ {x)ijj^ {y)ijj'' (z) , where e*-''^ is the complete antisymmetric sym- 
bol, and the color indices i, j, k can take the values 1,2, or 3. Thus, the con- 
struction below applies only for Nc — 3. (The generalization to arbitrary Nc 

^ A similar approach was adopted for the proton-proton scattering in Refs. [47,48]. 




0{x,y) 




d-'>%al-al){ai-al){al-al). 



(4.11) 
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is in principle possible, but the analysis becomes more complicated because 
the baryonic operator is then built with N^. quark fields.) By using the same 
technique as for dipole-CGC scattering, one obtains the following ^-matrix: 



Sr{x, y,z) = l^ e^"'e^^^{Vi{x)V^y)Vi{z))r, (4.12) 

where x, y, and z are the transverse positions of the three quarks. This op- 
erator is symmetric under any permutations of the three coordinates, and is 
normalized as St-{x, x, x) = 1. By using Eq. (4.4), it is easy to check that the 
odderon contribution is given again by the imaginary part of the S'-matrix : 

{B{x, y, z))r = Sr{x, y, z), (4.13) 
where the "3-quark odderon operator" B{x, y, z) has been defined as 



B{x,y,z) = 




(4.14) 



This is totally symmetric too, and satisfies the boundary condition B{x, x, x) — 
0, which is an immediate consequence of the normalization St{x, x, x) — 1. 

The simplest way to see that the operator (4. 14) is gauge invariant is to notice 
that this can be rewritten in terms of the manifestly gauge invariant operators 
shown in Eq. (2.14). Indeed, by using the identity^ 



y'^''e'^^Va{w)Vjm{w)Vkn{w) = detV{w) = 1, (4.15) 

where w is an arbitrary transverse coordinate, one can equivalently rewrite 
the 3-quark odderon operator B as 



+ triVj^V^V^V^V^V^) + triV^V^V^V^V^V^) - c.c. ] . (4.16) 

Note that this expression involves not only dipolar operators, but also higher 
multi-polar ones (quadrupoles and sextupoles). By construction, this expres- 
sion is independent off w when Nc — 3. Thus, it can be simplified by choosing 



^ We have already used a similar relation when we have specified the normalization 
of Sr{x,y,z). 
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w to be one of the quark coordinates, say w — z. Then, Eq. (4.16) reduces 
to: 




(4.17) 



which looks indeed considerably simpler, but where the symmetry in x, y and 
z is not manifest anymore (although, for Nf. — 3, we know that this expression 
Z5 totally symmetric, by construction). 

In particular, when two of the coordinates arc the same, the 3-quark odderon 
operator reduces to the dipole odderon operator, Eq. (4.3): 

B{x,z,z) = 0{x,z) = -B{x,x,z). (^"c = 3) (4.18) 

This is physically reasonable, because the diquark state is equivalent to an 
antiquark as far as color degrees of freedom are concerned, and can be most 
easily checked by setting y = z va. the r.h.s. of Eq. (4.17). Note however that, 
if one rather sets y — x va the same expression, then it is not immediately 
obvious that the ensuing expression for B{x, x, z) is indeed equal to —0{x, z) 
when A^c = 3, as it should. This is so because of the lack of manifest symmetry 
of Eq. (4.17), as alluded to above. Still, by diagonalizing the unitary matrix 
VjVz, and after performing some simple algebraic manipulations (relying on 
the fact that the corresponding eigenvalues Aj, i — 1,2,3, are pure phases 
and obey A1A2A3 = 1), it is possible to check that the equality B{x,x,z) — 
—0{x, z) — —B{x, z, z) holds indeed. 

In the weak field approximation, as obtained after expanding to lowest non- 
trivial order (i.e., to cubic order in a) the Wilson lines in any of the previous 
expressions for S(cc, t/, z), one finds again a gauge invariant linear combination 
of trilinear field operators with the color indices contracted with the d-symbol: 

(4. 

0} {« - al) + (at - tt^)} {{< - «y) + - < 

In fact, in this weak-field regime, the 3-quark C-odd operator is fully de- 
termined by gauge symmetry together with the requirement of total symme- 
try with respect to the external coordinates: Indeed, it can be checked that 
Eq. (4.19) uniquely emerges from the C-odd Green's function (5.1) after sym- 
metrization and imposing the finiteness condition (2.12). 



B{x, y, z) ~ 
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5 Odderon evolution in the dipole-CGC scattering 

In this and the next sections, we shall apply the general JIMWLK equation (in 
its dipolar form, cf. Eq. (2.11)) to the operators describing odderon exchanges 
constructed in the previous section, in order to deduce the evolution equations 
for the respective, C-odd, scattering amplitudes. 

We start with the simpler case of the dipole scattering, for which we shall 
discuss separately the weak-field limit (corresponding to a single scattering), 
and the general, non-linear, case (which includes multiple scattering). In fact, 
in this particular case, it is rather straightforward to write down directly the 
non-linear equations (see below) , from which the equations for the weak-field 
limit can be then simply deduced by linearization. Still, the more detailed 
approach that we shall follow below is instructive as a preparation for the 
more tedious case of the 3-quark operator, to be discussed in Sect. 6. 

5.1 Linear evolution and the odderon Green's function 

In the weak-field regime, where we can limit ourselves to a single odderon 
exchange, it is convenient to proceed as in Sect. 3 and introduce the (totally 
symmetric) odderon Green's function 

fr{x,y,z) = d'''^alalal)r, (5.1) 

in terms of which the weak-field version of the dipole odderon operator (4.10) 
can be rewritten in a form similar to Eq. (3.7) for the pomeron: 

0{x,y) ~ :^^{s(^frix,y,y) - fr{x,x,y)) + fr{x,x,x) - fr{y,y,y)}- 

(5.2) 

The discussion of the pomeron Green's function (3.8) in Sect. 3 appfies to the 
odderon Green's function (5.1) as well: The latter is not a gauge-invariant 
quantity, so its evolution under the original JIMWLK Hamiltonian would be 
afflicted by infrared singularities, which can however be regulated by using the 
(weak field version of the) dipolar Hamiltonian H^p. This yields a mathemat- 
ically well defined equation for fr{x, y, z), that we shall use as the definition 
of the odderon Green's function. 

Still as in the pomeron case, the equation obeyed by the scattering amplitude 
(5.2) — which is gauge-invariant — is not sensitive to the ambiguities which 
affect the definition of the Green's function, and comes out the same whatever 
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form of the Hamiltonian is used in its derivation. In fact, the only reason for 
introducing the Green's function (5.1) is the fact that it is for this function 
that we shall verify the BKP equation later on. 

Specifically, by using Eq. (3.3) for O — d'^'^a'^a^ya'^,, and after some lengthy 
algebra, one obtains 

—fr{x,y,z) = — d w- ^ ^ 

{frix, W, Z) + fr{w, y, z) - fr{x, y, z) - fr{w, W, z)) 

+ |2 cychc permutations I , (5.3) 

which is like applying the equation (3.9) for the 2-point Green's function to 
each pair of points within /^-(cc, z). It can be checked as before that this 
equation is well defined both in the infrared and in the ultraviolet. It is now 
straightforward to construct the equation obeyed by the linear combination 
in Eq. (5.2), and thus find that this is precisely the BFKL equation (3.10) 



X ((0(a;, z))r + (0(z, y)), - (0(x, y)),) (5.4) 

(we have also used the identity feagfdbgdabc = ^dedc), in agreement with an 
original observation by Kovchegov, Szymanowski and Wallon [43]. In view 
of the formal difference between the pomeron and the odderon operators, as 
given by Eqs. (3.7) and, respectively, (4.10), it may appear as a surprise that 
they obey both the same evolution equation. But this becomes more natural 
if one remembers that N and O arc, respectively, the imaginary part and the 
real part of the scattering amplitude T (with S = 1 + iT), and that in the weak 
field regime T obeys a linear equation with real coefficients — the linearized 
version of the first Balitsky equation, Eq. (2.9) — which is therefore separately 
satisfied by its real and imaginary parts. 

But since the initial conditions corresponding to N and O are different (in 
particular, they have different C-parities), so are also the respective solutions, 
and their behaviors at high energy. In Appendix A, we compute within the 
CGG formalism the C-odd initial conditions for some simple targets: a bare 
quark and a qq dipole. Namely, if the target is a single quark with transverse 
position Xq, one obtains 
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whereas for the more interesting case of a dipolar target (with the quark being 
at Xq and the antiquark at yo), one rather finds 



(0(..«)>..„^g '^-y-^ in3 ;: (5.6) 



c 



\X — Xq\ 


\y - 


yo\ 


\x - yo\ 


\y- 


Xo\ 



These expressions are in agreement with the corresponding results in Ref. [43] 
up to an overall numerical factor. As expected, these initial conditions are 
antisymmetric under the exchange of x and y, and thus satisfy the boundary 
condition (4.9). It is easily checked that this boundary condition is preserved 
by the evolution according to Eq. (5.4), as necessary for this equation to be 
well defined. 

The high-energy behavior of the odderon solution to Eq. (5.4) has been an- 
alyzed too in Ref. [43], where it has been shown that the projection of the 
general BFKL solution onto C-odd initial conditions selects (C-odd) BFKL 
eigenfunctions whose maximal intercept is equal to one. It turns out that these 
are the same eigenfunctions that were previously identified, by Bartels, Lipa- 
tov, and Vacca [12], as exact solutions to the momentum-space BKP equation 
[4,5]. Thus, in contrast to the pomeron solution to the BFKL equation, which 
at high energy rises exponentially with Y (N(Y) ~ e*^"''"^''^, with the BFKL 
intercept aw — I + {4:\n2)asNc/7r [1]), the corresponding odderon solution 
rises only slowly, as a power of y ~ In s. 

But, of course, these types of high-energy behavior (for either or O), which 
are mathematical consequences of the BFKL equation, are physically accept- 
able only so far as this equation is a correct approximation, that is, within 
the limited range of energies where the unitarity corrections are indeed neg- 
ligible. For higher energies, the evolution is governed by more complicated, 
non-linear, equations, to which we now turn. 



5.2 Non-Linear evolution 



For dipole-CGC scattering, the general evolution equations obeyed by the 
average amplitudes {N{x,y))r and {0{x,y))T in the strong field regime can 
be easily inferred from the first Balitsky equation (2.9) : Since the operators 
N{x, y) and 0{x. y) are, respectively, the real part and the imaginary part of 
the dipole 5-matrix S{x,y) = (l/iVe)tr(Kj\/y), cf. Eqs. (4.3) and (4.7), it is 
clear that the respective equations can be simply obtained by separating the 
real part and the imaginary part in Eq. (2.9). One thus obtain: 
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d_ 



[X - y) 



27r J {x — zfi^z — y) 

x{0{x,z) + 0{z,y)-0{x,y) 

- 0{x, z)N{z, y) - N{x, z)0{z, y))^, 

2nJ {x-zY{z-yY 

X (n{x, z) + N{z, y) - N{x, y) 

-N{x,z)N{z,y)+ 0{x,z)0{z,y) 



(5.7) 



(5.8) 



As is generally the case for the Balitsky equations, the equations above do not 
close by themselves, but rather belong to an infinite hierarchy. Interestingly, 
the non-linear terms in these equations couple the evolution of C-odd and 
C-even operators. For instance, the last term, quadratic in O, in the r.h.s. of 
Eq. (5.8) for {N)t describes the merging of two oddcrons into one pomeron. 
This process has not been recognized in previous studies of the Balitsky hi- 
erarchy, but the vertex connecting one pomeron to two odderons has been 
already computed in lowest order perturbation theory [26], and it would be 
interesting to compare such previous results with the corresponding vertex in 
Eq. (5.8) (which is essentially the dipole kernel). 



But the odderon-pomeron coupling which turns out to have the most dramatic 
consequences is the one encoded in the last terms in Eq. (5.7) for {0)t : As 
we shall shortly argue, this coupling leads to a rather rapid suppression of the 
C-odd contributions to scattering in the high energy regime where unitarity 
corrections start to be important (i.e., where {N)t- ~ 0{1)). To construct the 
argument without having to resort to an infinite hierarchy of equations, we 
shall restrict ourselves to the mean field approximation in which the non- 
linear terms in Eqs. (5.7)-(5.8) arc assumed to factorize. In the strong field 
regime, which will be our main focus below, we expect this approximation to 
be qualitatively correct [40]. 



In this mean field approximation, Eqs. (5.7)-(5.8) reduce to a closed system 
of coupled, non-linear, equations for {N)t- and {0)r ■ 



d_ 



{Oix,y)), 



a, 



I 



(fz- 



[x-y] 



27r J " ' {x — zy{z — yY 
X [{0{x, z))r + {0{z, y))r - {0{x, y))r 
- {0{x, z))r{N{z, y))r - {N{x, z))r{0{z, y)). 



(5.9) 
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d_ 




(5.10) 



X {N{x, z))r + {N{z, y))r - {N{x, y)), 
- {N{x, z))r{N{z, y))r + {0{x, y)); . 



The first of tliese equations lias been already proposed in Ref. [43], as a 
plausible non-linear generalization of Eq. (5.4). As for Eq. (5.10), this is the 
Kovchegov equation [39] supplemented by a new term describing the merging 
of two odderons. 

The Kovchegov equation has been extensively studied over the last few years, 
both analytically and numerically, and although the exact solution is not 
known, its general properties are by now well understood [49,50,51,52,53,54,55,56,57,58,59] . 
In its most synthetic description, due to Munier and Peschansky [59] (sec also 
[50]), this solution can be viewed as a front connecting the saturation regime 
where ^ A*" = 1 — this is reached for dipole sizes r = [a; — t/| much larger 
than the saturation length 1/Qs{t) — to the unstable regime at r <^ 1/(5s(t), 
where the amplitude is small {N <^ 1), but it rises rapidly with r, according to 
BFKL equation. When increasing r, the front propagates towards smaller val- 
ues of r (or larger transverse momenta), and its instantaneous position defines 
the saturation momentum Qs{t). The latter is found to rise exponentially with 
r: Q1{t) = Qle'^""'^, with c a numerical constant determined by the BFKL 
dynamics [17,56,57,58]. 

As we shall argue below, this mean-field picture of the pomeron exchanges 
is not significantly modified by the odderon contribution to Eq. (5.10). In 
particular, the values N — 1 and N — remain as (stable and, respectively, 
unstable) fixed points of the evolution described by Eqs. (5.9)-(5.10), but to 
them one should add a new fixed point, namely = 0, which is asymptotically 
approached at high energy. 

It is first easy to check that A^ = 1 and = are indeed fixed points at high 
energy. To also see that this is the only combination of fixed points in this 
limit, notice from Eq. (5.9) that, when increasing A^, the non-linear terms in 
this equation act towards suppressing the odderon [43]. It is thus consistent 
to assume that, at high energy, the odderon contribution represents only a 
small perturbation to the Kovchegov equation, so that the pomeron saturates 
in the standard way: A^(r, r) ~ 1 for r S> 1/Qs{t). Also, the dominant energy 
behavior of the saturation momentum, i.e., the value of the saturation expo- 
nent c, should not change, since this is fully determined by the linear (BFKL) 
part of the equation for A^. Using this assumption, it is possible to study the 
approach of O towards zero, and thus check the consistency of our hypothesis. 

^ Until the end of this section, we shall mostly use the simplified notations A'^ and 
O for the respective expectation values. 
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Namely, for sufficiently high energy, such that r 3> 1/Qs{t), the integral in the 
r.h.s. of Eq. (5.9) is dominated by relatively large dipoles, for which iV ~ 1. 
In this regime, the non-linear terms in this equation precisely cancel the first 
two linear terms there, and the equation simplifies to 



n C ri 7 

— {0{x,y))r^-as J ^{0{x,y))r = -as\n[Ql{Ty]{0{x,y))r,{5.11) 
which together with ln[(5^ (r)r^] = cas{T — Tq) immediately implies: 

{0{x,y))rc^exp!^~al{r-ro)^^{0{x,y))r, for r » l/g,(T). (5.12) 

As anticipated, this is a rapidly decreasing function of r, which is actually 
the same as the function describing the approach of the real part of the S- 
matrix towards the 'black-disk' limit S = (the Levin- Tuchin law) [49,37]. 
We expect the fluctuations neglected in the mean fleld approximation leading 
to Eqs. (5.9)-(5.10) to modify (actually, decrease) the value of the overall 
coefficient in the exponent of Eq. (5.12), but preserve the above qualitative 
picture [40]. 



6 Odderon evolution in the scattering of the 3— quark system 

The new feature which makes the 3-quark system conceptually interesting is 
the fact that the corresponding scattering amplitude depends upon three inde- 
pendent transverse coordinates. Therefore, already the lowest-order odderon 
amplitude, as given in Eq. (4.19), involves configurations in which the three 
exchanged gluons arc attached to different quark legs, and which are thus 
probing the complete functional dependence of the odderon Green's function 
defined in Eq. (5.1). As we shall further argue in Sect. 7, this in turns implies 
that the 3-quark — CGC scattering is a good theoretical laboratory to study 
the general solution to the BKP equation. 

Our discussion in this section will be restricted to the weak-field version of 
the 3-quark odderon operator, Eq. (4.19), which describes scattering via the 
exchange of a single odderon. This is indeed sufficient to discuss the cor- 
respondence with the BKP equation in the next section. Within the CGC 
formalism, there is no difficulty of principle (other than the tcdiousncss of the 
corresponding algebraic manipulations) which would prevent us from deriving 
the evolution equations satisfied by the general, non-linear, 3-quark opera- 
tors in Eqs. (4.14) or (4.16). However, these general equations are complicated 
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and not very ilhiminating: Through their non-hnear terms, they couple the 
3-quark operator to other color structures. This is especially manifest if one 
uses the form (4.16) of the 3-quark amplitude: this involves operators of var- 
ious multi-polar orders, which in the Balitsky hierarchy are coupled to other 
operators of even higher multi-polar moments. 

By using the expression of B{x, y, z) as a linear combination of odderon 
Green's functions, as manifest in Eq. (4.19), together with the equation (5.3) 
obeyed by the latter, one can deduce after a straightforward but lengthy cal- 
culation the following linear evolution equation for {Bxyz)T = {B{x, y, z))r' 



^ {j\^xwz) T ~l~ {-^wyz) r {-^xyz^r 



-{B 

+ (2 cyclic permutations). (6.1) 

(Of course, the same equation could be obtained also directly from Eq. (3.3) 
with O = Bxyz, but the corresponding manipulations would be even more te- 
dious; the introduction of the Green's function (5.1) at intermediate steps has 
also the advantage to better organize the calculation.) Note that Eq. (6.1) is a 
closed equation for {Bxyz)T, which was expected in view of gauge invariance: 
the only gauge invariant operators available are B^yz and Oxy — Bxyy (cf. 
Eq. (4.18)). 

Consider the structure of Eq. (6.1) in some detail: The linear combination 
of -B's in the integrand vanishes at the points w — x and w — y where lie 
the poles of the dipole kernel, so the latter are again innocuous. Also, one 
can easily check that the above equation is consistent with the relation (4.18) 
between the dipole and the 3-quark odderon amplitudes: if one sets z = y, 
Eq. (6.1) reduces indeed to Eq. (5.4) with Nc — 3. Lastly, Eq. (6.1) manifestly 
preserves the symmetry of the scattering amplitude under the permutation of 
its three coordinate variables. 

The initial condition can be calculated within the CGC formalism, in a similar 
way as for dipole-CGC scattering (cf . Appendix A) . For example, for a single 
quark at transverse position Xq, one finds 



{Bxyz') T=0 

(6.2) 

5 3 

As expected, this is symmetric under the permutation of the three external 
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coordinates, and satisfies the boundary condition {Bxxx)t — 0. 

In the next section, we shall argue that the high energy behavior of the solu- 
tion {B^y;^)r to Eq. (6.1) is controlled by the BLV solution [12] to the BKP 
equation, and thus has a maximal intercept equal to one. 



7 Comparison with previous approaches 

Within the traditional perturbative QCD approach to small-x evolution, the 
odderon exchange is viewed as the exchange of a composite object made of 
three reggeized gluons which evolves with energy according to the Bartels- 
Kwiecinski-Praszalowicz (BKP) equation [4,5]. From the point of view of the 
CGC formalism, this exchange is a single scattering, and should be compared 
to the weak-field limit of the corresponding CGC equations, as derived in 
the previous sections. In what follows, we shall demonstrate that the BKP 
odderon corresponds to the CGC Green's function introduced in Eq. (5.1): 
The BKP equation, which is traditionally written in momentum space, is 
essentially the Fourier transform of Eq. (5.1) for /r(a;,?/,z). Based on our 
previous discussion of gauge-invariant scattering amplitudes, we shall then 
argue that the natural Hilbert space for discussing odderon exchanges is in 
fact larger than the one which is generally used in the literature, and within 
which the BFKL Hamiltonian shows holomorphic separabihty [8,10]. 

7.1 Relation to the BKP equation 

The BKP equation is traditionally written in momentum space as 

^F,(fei, k2, fes) = ^ E / d'k[d'k',d% 6^'\k[ + k', + k',-q) 6^^\h - K) 

xHBFKL{ki-i, ki+i; k\^^)F^{k[, k'^, k'^), (7.1) 

where F-r{ki, k2, k^) is the Green's function for the exchange of three reggeized 
gluons (with transverse momenta fei, ^2, and k^, respectively) in a C-odd, 
color singlet, state. Furthermore, k^ = ki (and similarly k'^ = k[), and q is 
the momentum transfer. The factor | accounts for the fact that two of the 
gluons are in the color octet state. Finally, i^BFKL is the non-forward BFKL 
kernel including the virtual terms: 

iiBFKL(fci, k2, ki, k^) = klkl{ki - k[)^ 
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-n5^'\k[-k,)(ln^+ln^y (7.2) 



where e is an infrared cutoff which does not affect physical results. 

The function F^{ki, ^2, fcs) is often referred to as a 'scattering amplitude', but 
in fact the physical amplitudes are obtained only after convoluting this with 
appropriate impact factors. A generic odderon amplitude reads: 

0{q) = j (fk^(fk2d^k^ S^^Hki + k^ + h - q) 

X $proj(fci, k2, ks) Fr{ki, fca, fca), (7.3) 

where gauge symmetry requires the projectile impact factor $proj(fci, fe2, ^3) 
to vanish when ki — for some i (indeed, a zero-momentum gluon 'sees' 
the projectile as a whole, and the latter is globally colorless). Note that, when 
writing the scattering amplitude as in Eq. (7.3), the impact factor of the target 
is implicitly included in the definition of F^-. This is in line with the general 
approach in this paper, where the target is a generic 'color glass condensate', 
and the CGC Green's functions hke those in Eqs. (3.8) or (5.1) include the 
relevant information about the target impact factor. 

In the case where the projectile is a proton, the odderon impact factor is 
expected in the form (see, e.g., [48]) 

$proton(/ci,fe2,/e3) = / xd^yd^ ^,^t^^{x , y , z)\\2 e''^'^ - e^^^y - e'''^ 
x{2e'''^y - q'^^^ - e''=2'^)(2e^'=3"^ - e^'"^^ - e^*'^^), (7.4) 

where | ^proton (^c, 2/, ■Sj)^ is the proton light-cone wavef unction, with cc, t/, 
and z denoting the coordinates of the three valence quarks relative to their 
barycenter. (The condition x + y + z = Ois implicit in the definition of 
I ^proton P-) The exponential terms within the parentheses correspond to all 
the possible attachments of the three exchanged gluons to the quark lines in 
the proton. Note that $proton(fei = 0, fe2, ^a) = 0, etc., as expected. 

Furthermore, if the projectile is a virtual photon, the odderon couples to the 
dipole component of the photon Fock space, so that: 



*7*(fci,fc2,fe3) = J dzd^r\'if^*{z,r)\'^ (7.5) 
x(e*'°ii — e~''°^5)( e**"^^ — e~**'^^)(e*'°^^ — e"**"^^), 

where r is the dipole size, and ^'^* is the light-cone wavefunction describing 
the dissociation of the virtual photon, and can be computed in perturbation 



33 



theory with respect to ckem [60] . 



To make contact between this more traditional approach and our previous 
resuhs in this paper, it is convenient to introduce first the momentum space 
version of the odderon Green's function introduced in Eq. (5.1). We thus 
define: 



in terms of which the scattering amphtudes (0(a;, y))r and {B{x, y, z))r can 
be rewritten as 



{0{x, y))r = / (fkAdPk^fr{k^, ^2, fca) (7.7) 

^(gifeias _ gifciy)( gife2a5 _ gife2y) ( gifcsa: _ gifesj/) 



and, respectively. 



{B{x,y,z))r 

(fk^(fk2dFk^U{ki,k2M){'^^^'"' - ^^'^ - ^^''') 
x(2e**'^^ — 6*'°^^ — e**'^'")(2 6'*°^^ — e**'^'" — 6**°^^) 



(7i 



Note the similarity between the exponential factors in these equations and 
those in Eqs. (7.4) and (7.5) for the impact factors. This reflects the fact that 
the couplings between the exchanged gluons and the quarks (or antiquarks) in 
the projectile have been explicitly included in our definition of the scattering 
amplitudes. By inspection of the previous equations, it becomes clear that 
the CGC Green's function /^(fci, ^2, fcs) should correspond to the function 
Fr{ki, ^2, fcs) of the traditional BKP approach. To be able to compare these 
quantities, one also needs the equation satisfied by fr{ki, k2, k^), which is 
obtained after taking a Fourier transform in Eq. (5.3). A lengthy but straight- 
forward calculation shows that the resulting equation is the same as the BKP 
equation (7.1) up to terms proportional to delta-functions 5^^^ (fcj) (i = 1, 2, 3), 
which are however irrelevant for the calculation of the scattering amplitudes, 
since they do not contribute to the convolutions in Eqs. (7.3) or (7.7)-(7.8). 
This demonstrates the equivalence between the two formalisms, in so far as 
the single odderon exchanges are considered. 
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1.2 Comments on the Hilbert space for odderon solutions 

Following the remarkable discovery by Lipatov [8] that the BFKL Hamiltonian 
exhibits holomorphic separability when restricted to functions which belong 
to the Mobius space — by which we mean the functions /^{x.y.z,...) which 
vanish when any two coordinates coincide with each other (frix, x, z, . . .) = 0, 
etc.) — , much effort has been devoted towards finding solutions to the BKP 
equation (and, more generally, to its generalization which describes the ex- 
change of n reggeized gluons) within this particular Hilbert space of func- 
tions. This situation is mathematically appealing since the restriction of the 
BKP equation to a given holomorphic sector describes a dynamical system 
which has a sufficient number of hidden conserved charges (three in the case 
of the odderon) to be completely integrable^° [9,10], and which in fact can be 
identified as the XXX Heisenberg model of spin s — [10]. 

However, as recently reiterated in Rcf. [28], a careful inspection reveals that 
the full BKP Hamiltonian contains extra delta functions Uke 6^'^\x — y) in 
addition to the separable Hamiltonian given in [8,9,10]. One can discard these 
delta functions and restore holomorphic separability by working in the Mobius 
space. But this is only natural in the case of the BFKL pomeron, i.e., for the 
2-point Green's function /r(£C, J/), since in that case one can always redefine 



which ensures that /^(a;, x) — 0, without affecting the calculation of scattering 
amplitudes : The subtracted terms in Eq. (7.9), being independent of one of 
the coordinates, give vanishing contributions when convoluted with a colorless 
impact factor. However, for the odderon problem (n = 3), the restriction to 
the Mobius space is a highly nontrivial issue — it cannot be simply achieved 
via subtractions which preserve the physical amplitudes — , and therefore has 
imphcations on the physical relevance of the various solutions to the BKP 
equation. 

And, indeed, among the two exact solutions to the (odderon) BKP equations 
which do not vanish rapidly at high energies [11,12], only one of them — the 
Janik-Wosiek (JW) solution [11], which has an intercept slightly lower than 
one — lies indeed in the Mobius representation, e.g.. 



For the generalization of the BKP equation to a system of n reggeized gluons, 
with n > 3, the same property holds only in the large-iVc Umit in which one can 
restrict the BFKL-like pairwise interactions to neighboring gluons [8,9,10]. 

In fact, a brief comparison with Eq. (3.7) reveals that the subtracted 2-point 
Green's function can be identified with the pomeron scattering amplitude. 




1 



1 



(7.9) 
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f^'^{x,x,z)^0, 



(7.10) 



and has been constructed by exploiting conformal symmetry and integrabil- 
ity. But, clearly, this solution does not couple to a dipole in the present, 
leading-logarithmic, approximation. Besides, although in principle it can cou- 
ple to a 3-quark system, the corresponding amplitude constructed according 
to Eq. (4.19) (where we identify d"''^'^ (a'^aya^.) r = fr^ix, y, z)) has the rather 
curious property to vanish at equal points {{B{x, z, z))r = 0, etc.), for which 
there is no compelling physical justification. For instance, the initial condition 
(6.2) does not have this property. Thus, the solution to Eq. (6.1) corresponding 
to this particular initial condition will not belong to the Mobius space. 

Similarly, the other exact solution known for the odderon BKP equation, due 
to Bartels, Lipatov, and Vacca (BLV) [12], does not belong to the Mobius 
space either. For this solution, the dominant intercept is exactly one, so the 
corresponding amplitude has only a weak energy dependence (weaker than any 
power of the energy). As explicitly verified in Ref. [43], using the BLV solution 
within Eq. (7.7) for the dipole odderon amplitude (0(cc,t/))^ is equivalent to 
constructing the general C-odd solution to the BFKL equation (5.4). Thus, 
the BLV solution appears as the physical Green's function describing the odd- 
eron exchange between a virtual photon and some generic target (in the dilute 
regime where saturation effects in the target are unimportant). 

Furthermore, by inserting the BLV solution into Eq. (7.8), one obtains a par- 
ticular solution to Eq. (6.1) for the amplitude {B{x, y, z))r describing odderon 
exchanges between a 3-quark system and a dilute target . The following ar- 
gument suggests that this particular solution should in fact yield the dominant 
behavior of (S)^ at high energy; that is, without actually solving Eq. (6.1), 
one can conclude that the intercept for {B)^ is exactly one, and is determined 
by the BLV solution. The argument relies on the relation (4.18) between the 
odderon amplitudes for the 3-quark system and that for the dipole. Let l + u) 
denote the intercept for (-B)r, that is. 



where h is only slowly varying with r (slower than any exponential). Setting 
y — z, one finds 



Since, on the other hand, (0(cc, z))^ must be a C-odd solution to the BFKL 
equation (5.4), the analysis in Ref. [43] demonstrates that the largest possible 

To fully determine the physical amplitude, one must still convolute the BLV 
Green's function with the target impact factor. 



{B{x,y,z))r = h{x,y,z;T)e' 



(7.11) 



{0{x,z))r = h{x,z,z;T)e' 



(7.12) 
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value of a; is a; = 0. This argument misses the component of the BKP solution 
which vanishes at equal points. However, the largest intercept of this kind, 
corresponding to the JW solution, is known to be less than one [11]. 



8 Conclusions 

In this paper we have given the first discussion of the odderon problem within 
the framework of the effective theory for the color glass condensate. The gen- 
uinely CGC part of the analysis has been rather straightforward: The operators 
describing C-odd exchanges between the CGC and two simple projectiles — a 
qq color dipole and a 3-quark system — have been constructed in the eikonal 
approximation, which includes multiple scattering to all orders via Wilson 
lines. By acting with the JIMWLK Hamiltonian on these operators, we have 
then deduced the evolution equations satisfied by the corresponding scattering 
amplitudes. This is a standard, but generally quite lengthy, procedure, that we 
have considerably simplified by introducing the dipolar form of the JIMWLK 
Hamiltonian. In particular, for a dipole projectile, we have recovered the equa- 
tions previously obtained in Ref. [43], that we have generalized here beyond 
the mean field approximation. 

What turned out to be more subtle, however, was the correspondence with the 
traditional perturbative QCD approach (in particular, with the BKP equation) 
in the limit where the scattering is weak. Although our equations become linear 
in this limit, they still apply to gauge-invariant scattering amplitudes, and not 
directly to Green's functions. Besides, they are a priori written in coordinate 
space. But the use of the dipolar version of the JIMWLK Hamiltonian has 
allowed us to write down a well-defined equation for the odderon Green's 
function, which after Fourier transformation to momentum space turned out 
to be the same as the BKP equation. 

Our analysis has emphasized the subtlety involved in the Fourier transforma- 
tion of the BKP equation, and the importance of the structure of the external 
probe for selecting physical solutions to this equation. While the latter point 
was already stressed in Refs. [12,43] in the context of the dipole scattering, 
the ability of our formalism to deal with a 3-quark system (which has three 
independent coordinates) has made this point even clearer, thus shedding light 
on the Hilbert space to be used in relation with the BKP equation. From the 
viewpoint of Eq. (6.1), the solutions which belong to the Mobius represen- 
tation (i.e., which vanish at equal points: {B{x,z,z))t = 0, etc.) are very 
special ones, and are unlikely to be realized during the evolution from phys- 
ical initial conditions. Rather, the solutions to Eq. (6.1), and therefore also 
to the BKP equation, which are singled out by our gauge-invariant 3-quark 
amplitude and the respective initial conditions do not vanish at equal points. 
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This observation, together with the relation Eq. (4.18), has led us to conclude 
that the highest odderon intercept for the 3-quark-CGC scattering is exactly 
one, so like for the dipole-CGC scattering [43], and is described again by the 
Bartels-Lipatov-Vacca solution [12] to the BKP equation. 

As mentioned in the Introduction, our weak field analysis of the odderon 
problem (n = 3, C = — 1) is intended as a first step in a systematic study of 
multireggeon (n > 3) exchanges in the CGC formalism, with the aim of clar- 
ifying the relation between this formalism and more traditional approaches, 
hke GLLA. Although successful in estabhshing the correspondence with the 
BKP equation, our previous analysis has also revealed a few subtleties which 
may lead to difficulties when trying to extend this approach to processes with 
more than three reggeons. We have seen indeed that the Green's functions in 
the CGC formalism make sense only as building blocks (in the sense of linear 
combinations) for scattering amplitudes in the weak field limit. Thus, in order 
to have a meaningful definition for a n-point Green's function (corresponding 
to a n-reggeon exchange), one must first identify appropriate gauge-invariant 
amplitudes whose weak-field expansion starts at order n in a"'{x). The poten- 
tial difficulty with this approach, however, is that there is a priori no guarantee 
that the nonlocality in x" inherent in the Wilson lines will disappear in the 
expansion leading to the Green's functions. Recall, for instance, our previous 
construction of the weak-field odderon operators, Eqs. (4.10) and (4.19), from 
the corresponding non-linear operators, Eqs. (4.3) and (4.14), respectively: In 
that case, the nonlocality in x" has disappeared from the final results only 
'accidentally', because of the presence of the totally symmetric tensor d"^'^. An 
alternative procedure which looks promising would be to construct the gauge- 
invariant linear combination directly in the weak-field limit, starting with a 
non-invariant Green's function and imposing on it the finiteness condition 
(2.12). It remains as an interesting open question whether any of the methods 
mentioned above can be used to construct arbitrary n-reggeon exchanges. We 
leave this and related issues for future work. 
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A C-odd initial conditions in the dipole-CGC scattering 



The gauge field a.% is created by a color source in the target CGC. It is 
given by (in the covariant gauge) 



where is an infrared cutoff which will disappear from the final results. For 
a single quark at transverse position Xq, the color source is given by 



p\z) = Q'^5'^^\z-Xo), (A.2) 
and for a dipole made of a quark at Xq and an antiquark at 1/0, 



= [5^^\z - xo) - S^'\z - yo)] , (A.3) 

where is the color charge Q"' = g J ipH'^ip. We obtain the initial conditions 
for the evolution equation by substituting the gauge fields created by these 
"unevolved" targets into Eq. (4.10) and (4.11), and taking the average over the 
"random" gauge field. For the unevolved targets the average over the random 
configuration simply reduces to the color average. More precisely, in order to 
evaluate the average (Q'^Q^Q^), we replace the c-number charge by a color 
matrix gf^ and take the trace -^ti {gf^ gt'^ gt^). Thus, it can be evaluated as 
follows: 



d<^bc(^QaQbQc^ — > d^^'^trigt" gt'' gf) ^ ^d'''"'hr{{e,t''}t'') 

3 3 
= A_rf«^c^«^c ^ _A_^^2 _ 4)^^2 _ ^A.4) 

This yields the results (5.5) and (5.6). The technique adopted here for the 
average was first proposed by lancu and Mueller for the onium-onium scat- 
tering to show the equivalence between the color dipole picture and the CGC 
formalism [45]. In fact, this procedure goes beyond the original McLerran- 
Venugopalan (MV) model. This is because the MV model is formulated with 
the Gaussian random source and the average of odd number of sources such as 
{p°-p^p'^) is simply vanishing. In order to describe the odderon, it is necessary 
to extend the MV model so that there is nontrivial correlation among three 
sources. On the other hand, such nontrivial correlation is correctly encoded 
in the CGC framework, as we claim in the present paper. Indeed, the dipole 
JIMWLK equation (3.3) in the weak-field limit suggests a similar equation 
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for the weight function iyr[«], but the Gaussian weight function is not the 
exact solution to this evolution equation. It arises only in the mean-field like 
approximation discussed in Ref. [37]. 
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